Exact Charge Dynamics in the Supersymmetric t-J Model 
with Inverse-Square Interaction 
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Dynamical charge structure factor N(Q, u>) with Q smaller than the Fermi wave number is derived 
analytically for the one-dimensional supersymmetric t-J model with 1/r 2 interaction. Strong spin- 
charge separation in dynamics is found. N(Q, ui) with given electron density does not depend on the 
spin polarization for arbitrary size of the system. In the thermodynamic limit, only two holons and 
one antiholon contribute to N(Q,uj). These results together with generalized ones for the SU^T) 
supersymmetric t-J model are derived via solution of the Sutherland model in the strong coupling 
limit. 
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Strong correlation between one-dimensional electrons 
has drastic consequences to their physical properties. 
The most spectacular is the spin-charge separation; the 
elementary excitations are spinons and holons both of 
which obey the fractional statistics [Q. This is in con- 
trast with the fermionic Landau quasi-particles realized 
in higher dimensions. In the supersymmetric t-J model 
with 1/r 2 interaction , these fractional particles appear 
in the simplest manner. Exact thermodynamics for the 
model BQl can be interpreted in terms of free spinons 
and holons. The charge susceptibility, for example, de- 
pends on the electron density, but is independent of the 
spin polarization. This kind of independence is referred 
to as the strong spin-charge separation Hlj. It is nat- 
ural to ask how the spin-charge separation appears in 
dynamical correlation functions. 

In the absence of charge degrees of freedom, the super- 
symmetric t-J model is reduced to the spin chain called 
the Haldane-Shastry model ||||. For the latter model 
Haldane and Zirnbauer j?j derived exactly the dynami- 
cal structure factor which is non-zero only over a com- 
pact region in the momentum-frequency plane. Such re- 
gion is called the support. For the supersymmetric t-J 
model, Ha and Haldane (^] analyzed numerical results for 
finite-size systems, and found that only a few number of 
elementary excitations contribute to spectral functions, 
namely dynamical structure factors. They proposed a 
compact support for each spectral function in the ther- 
modynamic limit, but did not obtain the spectral func- 
tions themselves. Analytical knowledge should obviously 
provide deeper insight into strong correlation effects on 
dynamics. Recently, Kato partially derived a spectral 
function corresponding to the removal of a particle from 
the high-density limit So far all analytical work in 
the literature does not include charge degrees of freedom 
in the ground state, and derivation of charge dynamics 
has remained a theoretical challenge. 

In this Letter, we report on exact analytical results on 



charge dynamics for the one-dimensional supersymmetric 
t- J model with 1/r 2 interaction at zero temperature, and 
its generalization from SU(2,1) to SU(-ftT,l) supersymme- 
try. We show explicitly how the strong spin-charge sep- 
aration manifests in charge dynamics. The charge cor- 
relation function is relevant to dielectric and elastic re- 
sponses of real systems with quasi-one-dimensional elec- 
trons. Our work should motivate further analytical work 
toward other dynamical correlation functions which are 
relevant to various spectroscopies such as neutron scat- 
tering and photoemission. 

We consider the supersymmetric t-J model given by 
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where ci G is the annihilation operator of an electron at 
site i with spin a, rii is the number operator and Si is 
the spin operator. The projection operator V excludes 
double occupation at each site. The transfer energy tij 
and exchange one Jij are given by tij = Jij/2 = tD~^ 
where D L j = (N/tt) sin (it (i — j) /N) with N being the 
number of lattice sites. This model is generalized to the 
SL^if,!) supersymmetry as 
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Here we have introduced a graded permutation operator 
Pij = J2 a (3 X'* 13 Xj a 0p, where X^ a changes a state a 
to (3 at site j, with a, (3 being either (hole state), or 
one of a = 1,- • - ,K (spin state). The sign factor Op is 
— 1 if j3 = and 1 otherwise. In the case of K = 2 the 
expression (||) is reduced to eq.(||). 

In terms of the density fluctuation operator uq = 

N ~ 1/2 Ef=i Ezli x r exp(iQl) with Q > 0, the dynam- 
ical charge structure factor is defined by 
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N (Q, u) = IH™Ql°>| 2(5 {u;-E u + E q ) 
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where \v) denotes a normalized eigenstate of the Hamil- 
tonian (Q) with energy E v . 

In the SU(if,l) t- J model, a holon has the spectrum 
e h(q) = Ktq[2n(l — n)/K + q}/2 with n being the average 
electron number per site, and another charge excitation 
called the antiholon has the spectrum e a (q) = £g[27r(l — 
n) — q]/2 |||§. The latter represents the excitation from 
the condensate of holons. We first give our main result 
of this Letter, and then its derivation. The system has 
N a electrons with spin er. In the small momentum re- 
gion < Q < fc£? in , with fc^? in being the smallest Fermi 
number of electrons, i.e., kf in = 7rMin(iVi, • • • N K )/N, 
we derive N(Q,u>) in the form: 
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where A„ = K/(2%) ]jf =1 T(l/K)/T(j/K) 2 . Surpris- 
ingly this N(Q ) u>) has almost the same form as that of 
the spinlcss Sutherland model with the coupling param- 
eter 1/K |lOUllf| , In the latter model excitations consist 
of K particles created above the pseudo Fermi level with 
charge — e, and a hole with charge Ke. The only differ- 
ence in N(Q,oj) between the two models is the integra- 
tion range of momenta. Thus the result above can be 
interpreted in terms of the fractional statistics JlJ] . 

With K = 2, it can be shown that N(Q,u>) diverges 
as [eh(Q) — lo]^ 1 / 2 as the frequency approaches the up- 
per edge corresponding to the holon dispersion ||l3[| . The 
static structure factor N(Q) — f du>N{Q,iS) reproduces 
the known result juj obtained for the Gutzwiller wave 
function. This should be so since the exact ground state 
of the t- J model is given by the Gutzwiller wave function 

As in the previous study of thermodynamics 
|p|,^| JT5|Jl^| , we use the Sutherland model with the 
SU(ii',l) internal symmetry as an auxiliary, and take the 
limit (3 — > oo of the coupling parameter [3. Then the 
particles crystallize with equal distance from each other, 
and the resultant dynamics excluding phonons and uni- 
form motion of the center of gravity is mapped to that 
of the SU(-?Sf,l) t-J model given by eq.(^|). The merit of 
using the Sutherland model is that much more is known 
about properties of the eigenfunctions than those for the 
t-J model. The Sutherland model with arbitrary internal 
symmetry is written in the form 
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where L is the size of the system, and Sij represents the 
exchange operator of coordinates of particles i and j [ [l7| . 
In order to reproduce the lattice model we take the limit 
of large (3 and M, keeping the ratio t = (3/M fixed. Here 
we take the lattice parameter L/N as the unit of length. 
Note that the symmetry of the wave function leads to the 
relation = —Pij within the Hilbert space of identical 
particles. 

Let us describe derivation of N(Q,u>) taking the 
SU(2,1) Sutherland model. Generalization to the 
SU(if,l) supersymmetry is straightforward. The system 
has iVh holes, iVj up-spin electrons and N± down-spin 
ones, whose coordinates are represented by x\ for i-th 
hole, x\ for i-th up-spin electron and x\ for z-th down- 
spin electron. We arrange them as x = (xi, X2, ■ ■ ■ , xn) — 
(x±, • • • , x\ h , x\, ■ ■ ■ , xjy^ , x~i , • • • , xjy^ ) = (x^ 1 , x^ , x^) and 
assume that N b is even, and both Nf and iVj_ are odd. 
Then the ground state is non-degenerate. The interval 
/ = [1,N] denotes {i 6 Z|l < i < N} for a positive in- 
teger N. We define I h = [1, N b ], J T = [N h + l,N h + iV T ] 
and 1 1 = [Nh + Nf + 1, N] . Without loss of generality we 
assume N± > Np The wave function of the ground state 
for a set of N b , N^,Ni is given by, 
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where the complex coordinates z = (z\, ■ ■ ■ , zjv) are re- 
lated to the original ones x — (xi,---,xn) by Zj — 
exp(2itixj/L). 

The spectrum of the Sutherland model is conveniently 
analyzed with use of a similarity transformation gener- 
ated by O = Y[ l<jeI (1 - Zjj Zi f zf l ~ 1] ' 2 - The 

transformed Hamiltonian H = 0~ 1 HO — eo describes the 
excited states relative to the spinless particle model with 
the ground-state energy e = n 2 /3 2 N(N 2 -l)/6ML 2 . We 
obtain 
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Here di is called the Cherednik-Dunkl operator Eqj and 
is given by 

di = Zi- - — — (1 - Sik) 
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It is known that di can be diagonalized simultaneously 
with real eigenvalue A^. The resultant eigenfunction are 
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polynomials E\(z; (3) of Zi such as Yl i + . . ., and are 
called non-symmetric Jack polynomials Jl§| ]. Since we 
are dealing with identical particles, eigenfunction should 
satisfy the following conditions of the SU(2,1) supersym- 
metry: 

(i) symmetric with respect to exchange between z^'s; 

(ii) anti-symmetric with respect to exchange between zf's 
with the same a. 

By taking linear combination of E\(z;(3), we can 
construct a polynomial K\{z;8) with the SU(2,1) su- 
persymmetry. The latter belongs to a family called 
the Jack polynomial with prescribed symmetry [jl9| pl[ . 
We specify the set of momenta as (A^A^A^) = 
(A^---,A^ h) Al ) --- ) A^ T ,Ai,---,Ai ri ) with Ai > ••• > 

^AV^l > ■•■ > ( a = T;i)- ln tms wa y we can 

parameterize K\{z\8) by using A = (A h ,A T ,A^). At 

, /x^) with 
.(JVi - l)/2), p) = 
■■,(N l -N 1 )/2) and 



the ground state we have A = Acs = (p h ,P' 
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We define the inner product for functions f(z) and g(z) 
in iV-complex variable z = (zi, ■ ■ ■ , zn) as follows: 
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where f{z) denotes the complex conjugation of f{z). 
As is clear from the definition of the transformation O, 
K x (z;8) is orthogonal with respect to the above inner 
product. In order to derive the norm of K\(z;fi), we 
generalize the procedure taken in rcfs. |f2l|-p3||. The re- 
sult of lengthy calculation is given by 



(K x ,K x ) =N h \N l \N^ 1 (E x ,E x ) , 
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where (E x , E x )q denotes the norm of the non-symmetric 
Jack polynomials. We refer to the literature |L8| for the 
explicit form of the norm since it requires many lengthy 
combinatorial quantities. The quantity p x = p x P\P\ is 
given by 
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In the lattice model, the density correlation of parti- 
cles is the same as that of holes because of the complete- 
ness relation X) ff=T j_ Xf 7 + X i° = L In the Suther- 
land model, on the other hand, the correlation function 
of the sum of hole-density and particle-density gives the 
phonon spectrum. It can be shown that the intensity of 
the phonon correlation function is smaller than the spin 
correlation function by a factor of 0(/3 _1 ). In the lead- 
ing order in 8~ 1 , the correlation function of holes is the 
same as that of particles. In view of this situation we 
consider the power sum p m = Yliei Z T f° r holes. The 



hole-density hq can be expressed as hq — p m /\/N with 
Q = 2nmjN ]l0| , p| . In order to calculate N(Q,u>), we 
need to know the expansion coefficient c\ which appears 
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It is difficult to derive c\ for general values of m. How- 
ever, there occurs drastic simplification for the small mo- 
mentum region with m < (iVj — l)/2. In this region one 
can show that Pm^ gs has only charge excitations due to 
the triangular structure of the polynomial K\(z; 0) p3. 

Without spin excitations, we have the explicit expres- 
sion of K\(z;P) following Baker and Forrester [jl9| as 
follows: 
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where a — 2 + 1/ (3 and J"(zh) = J" (z±, ■ ■ ■ , z N=h ) is the 
symmetric Jack polynomial iQ. The charge excitations 
are characterized by the set v = (y\ , • • • , i>N h ) > which is 
related to A by A = (z/, 0, 0) + Acs- 

Since p m can be expanded by the symmetric Jack 
polynomials [p5|, we can determine the coefficient c\ in 
eq.(|l^). With use of a square s = in the Young 

diagram v, we introduce the notations 

[Or- II K^)-C(s)], h« = Hla(a u (s)+l)+l v (s)] 



where a„(s) = Vi — j is the arm- length, a'„(s) = j — 1 is 
the arm-colength, l v (s) = v'j—i with i/j the length of the 
column j is the leg- length, and l' u (s) = i — 1 is the leg- 
colength. We then obtain the dynamical structure factor 
of the Sutherland model with the coupling parameter (3 
as follows: 
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where Q = 2-Km/N and \v\ — J^f^i^j- The excitation 
energy is given by 

8 Nh 

Ex - E Xgs = Jj^ ]T[a^ 2 + (N h + 1 - 2i) Vi ]. (15) 

i—l 

By taking the limit 8 — > 00 with the ratio t = 8/M 
kept fixed, we can derive N{Q,u>) of the t-J model (Q). 
In this limit N(Q,u>) simplifies dramatically by the fol- 
lowing fact: The coefficient [0']" in eq.(pH) vanishes if 
the partition v contains s = (3,2). Then N(Q,ui) is 
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determined by three parameters which are related di- 
rectly with momentum of the elementary excitations: two 
holons and one antiholon. We emphasize that N(Q,u>) 
does not depend on the spin polarization as long as Q 
is smaller than the smallest Fermi wave number of all 
electrons. 

Derivation of dynamical charge structure factor for 
SU(K,1) model proceeds in a similar manner. We 
write the momenta of K holons and one antiholon by 
Ai, • • • , \k and A', respectively. Then N(Q,u) for finite 
iV and arbitrary K can be expressed as follows: 

N{Q,u) = B m 2 S mm S(u-E xx >) 
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where ||A|| = A' + £jLi(Aj - 1) and B = 
2ttAqK 1 ^ k Kl/N . The excitation energy is given by 

E X y= t -(^ 2 {X>(N h + K-X>) 

K 

+ £[(A* - l)(JV h + K(Xj + 1) - (2j - 1))]}. (17) 

3=1 

In the special case of Xk = where less than K holons 
are excited, we need to modify the form of eq.(16). How- 
ever, this case can be neglected in the thermodynamic 
limit. We have checked the validity of eq.(jl6|) together 
with the modified one in the case of K = 2 by comparison 
with numerical results up to N = 16. 

Following the same procedure as that in the spinless 
Sutherland model we obtain the expression (Q) in 

the thermodynamic limit. The present derivation makes 
it clear that the spin-charge separation persists in sys- 
tems with finite size. In the SU(if,l) model, a in eq.jl^) 
is given by a = K + /3 , which tends to K in the limit 
of (3 — * oo. Since 1/a can be regarded as the coupling 
constant in the charge sector of the fractional statistics, 



N(Q,u)) given by eq.(Q) has the apparent coupling con- 
stant 1/K. 

The authors would like to thank Y. Kato for valuable 
discussions. T.Y. and Y.S. wish to acknowledge the sup- 
port of the CREST from the Japan Science and Technol- 
ogy Corporation. T.Y. wishes to thank the support of 
the Visitor Program of the MPI-PKS. 



[1] F.D.M. Haldane, Phys. Rev. Lett. 66, 1529 (1991). 
[2] Y. Kuramoto and H. Yokoyama, Phys. Rev. Lett. 67, 
1338 (1991). 

[3] Y. Kuramoto and Y. Kato, J. Phys. Soc. Jpn. 64, 4518 

(1995) . 

[4] Y. Kato and Y. Kuramoto, J. Phys. Soc. Jpn. 65, 1622 

(1996) . 

[5] F.D.M. Haldane, Phys. Rev. Lett. 60, 635 (1988). 

[6] B.S. Shastry, Phys. Rev. Lett. 60, 639 (1988). 

[7] F.D.M. Haldane and MR. Zirnbauer, Phys. Rev. Lett. 

71, 4055 (1993). 
[8] Z.N.C. Ha and F.D.M. Haldane, Phys. Rev. Lett. 73, 

2887 (1994). 
[9] Y. Kato, Phys. Rev. Lett. 81, 5402 (1998). 
[10] Z.N.C. Ha, Phys. Rev. Lett. 73, 1574 (1995); 74, 620 

(1995) (E); Nucl. Phys. B 435, 604 (1995). 
[11] F. Lesage, V. Pasquier and D. Serban, Nucl. Phys. B 

435, 585 (1995). 
[12] Y. Kato, T. Yamamoto and M. Arikawa, J. Phys. Soc. 

Jpn. 66, 1954 (1997). 
[13] ER. Mucciolo, B.S. Shastry, BD. Simons and B.L. Alt- 

shuler, Phys. Rev. B 49, 15197 (1994). 
[14] F. Gebhard and D. Vollhardt, Phys. Rev. B 38, 6911 
(1988). 

[15] A. Polychronakos, Phys. Rev. Lett. 70, 2329 (1990). 
[16] B. Sutherland and B.S. Shastry, Phys. Rev. Lett. 71, 5 
(1993). 

[17] A. Polychronakos, Phys. Rev. Lett. 69, 703 (1992). 

[18] I.G. Macdonald, Symmetric Functions and Orthogonal 

Polynomials, (American Mathematical Society, 1997). 
[19] T.H. Baker and P.J. Forrester, Nucl. Phys. B 492, 682 

(1997) . 

[20] C.F. Dunkl, Comm. Math. Phys 197, 451 (1998). 
[21] Y. Kato and T. Yamamoto, J. Math. A: Math. Gen. 31, 
9171 (1998). 

[22] K. Takemura and D. Uglov, J. Phys. A: Math. Gen. 30, 

3685 (1997). 
[23] Y. Kato, private communication (1997). 
[24] Y. Kato and Y. Kuramoto, Phys. Rev. Lett. 74, 1222 

(1995). 

[25] P.J. Hanlon, R.P. Stanley and JR. Stembridge, Con- 
temp. Math. 138, 151 (1992). 



4 



